Spectral representation of some non 
stationary a-stable processes 



Nourddine Azzaoui 



Institut de Mathmatiques de Bourgogne (1MB) 
9 avenue Alain Savary 
21000 Dijon, France 



Abstract 



In this paper, we give a new covariation spectral representation of some non station- 
ary symmetric a-stable processes (SaS). This representation is based on a weaker 
covariation pseudo additivity condition which is more general than the condition of 
independence. This work can be seen as a generalization of the covariation spectral 
representation of processes expressed as stochastic integrals with respect to indepen- 
dent increments SaS processes (see Cambanis (1983)) or with respect to the general 
concept of independently scattered SaS measures (Samorodnitsky and Taqqu 1994). 
Relying on this result we investigate the non stationarity structure of some harmon- 
isable SaS processes especially those having periodic or almost-periodic covariation 
functions. 

Key words: Symmetric a-stable processes, Covariation spectral representation, 
Harmonisable processes, Periodically (almost-periodically) covariated processes. 



1 Introduction 

In this work, we are interested in a family of stochastic processes having infinite 
second order moments. It is about the class of symmetric a-stables (SaS) 
processes, with (1 < a < 2), having the stochastic integral representation: 
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where ^ is an SaS stochastic process. By using the concept of the covariation 
that will be recalled later in (7) or (8) and its relationship with the convergence 
in probability, Cambanis [l[ have given general conditions for the existence of 
the integral (1), they are: 

• (I) The process ^ is right continuous with respect to the convergence in 
probability. 

• (II) for all linear combination ( of increments of ^, the map v : t i — > 
[^(t): C]a is of bounded variations where [.,.]□, denote the covariation. 



Due to their impulsive nature, these processes provide appropriate models in 
various application fields: for example in financial and econometric modeling 
see Uchaikin and Zolotarev [2o|| , communications, signal processing see Nikias 



and Shao [17|] , nuclear physics, astronomy see Uchaikin and Zolotarev [20|] and 



references within... As an alternative to the covariance, which is undefined in 
this case, the covariation was introduced by Miller [l5|. Although this new 
dependence measure was conceived to substitute the covariance when 1 < 
a < 2, it is not as convenient because it does not have some of its suitable 
properties. One of its great defects is that it is not, in general, additive with 
respect to its second variable. This last property play an important role in 
integral representation of the covariation with respect to a spectral measure 
(control measure). In the case of symmetric a-stables processes, Cambanis 
have used the concept of independence as a sufficient condition for the 
covariation additivity what have enabled him to build a spectral representation 
of the covariation. This case have been widely studied in literature and was 
the keystone of several works in theoretical as applied fields. Among others. 



it was used in hnear regression (Samorodnitsky and Taqqu [l9|), prediction 
problems (Cambanis and Miamee [2|,[5|)and spectral analysis of stationary 
processes see Samorodnitsky and Taqqu 19|] and reference within. 



In this paper, we generalize works of Cambanis fl"] by replacing the incre- 
ments independence condition on the process ^, by a weaker one permitting 
the additivity of the covariation. Using this property we give, in the general 



sense of Morse and Transue [16| , a new integral representation of the covaria- 
tion function of the process X with respect to a characterizing bimeasure. This 
representation which is similar to the covariance spectral representation of sec- 
ond order processes, will play an important role in the study of a wide class 
of non- stationary a-stables processes. For instance, the class of second order 
harmonisable processes, introduced by Loeve, have been extensively studied 



by Rao [18|, Chang and Rao [6|. In particular, periodically and almost period- 
ically correlated second order processes have been thoroughly investigated, see 
for instance Hurd (lil : 13], Dehay [sj, Dehay and Hurd fq. In the Cambanis' 



s 

spectral representation, where ^ have independent increments, all harmonis- 
ables processes are covariation stationary: their covariation depend only on 



time difference, see Cambanis However, our spectral representation give 
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an understanding of the exact manner in which the increments of the processes 
^ are not independents gives some information about the non-stationarity na- 
ture of the process X. In particular, similarly to second order processes, we 
study an important class of harmonisables non stationary processes for which 
the covariation is periodic or almost periodic that we call {periodically and 
almost-periodically covariated processes). We show that some fundamental re- 
sults concerning second order periodically and almost periodically correlated 
processes remain true for our periodically and almost-periodically covariated 
a-stable processes. This paper is organized as follows: we begin by a brief 
overview on stable variables and processes and some preliminary results. We 
give, in theorem 2.1, our first result concerning a weaker condition for covari- 
ation additivity. After a brief recall on Morse and Transue integration, we 
provide our main spectral representation, finally, a classification of harmon- 
isable processes according to the structure of their characterizing bimesure is 
provided. All the proofs will be given in appendix. 



2 Definitions, notations and preliminary results 



A real centered random vector X'^={Xi, X2, . . . , X^) is symmetric a-stable 
(SaS) if and only if its characteristic function is given by: 



4>x4'ti, ...,td)= exp{- / 



d ° 

1=1 



C?r(Xi,...,Xd)(Sl, • • • 



(2) 



where T {^Xi,...,Xa) is an unique finite symmetrical measure defined on the unit 
sphere of M*^, see Samorodnitsky and Taqqu [l9|]. When X''- is a complex SaS 
vector, its characteristic function is defined, for all complex 6* = (6*1, . . . , 6*^) in 
by: 



<l>x(^i, . . . = E (exp{^7^e(^,X)}) = E (^exp{^7^e^^^fcXfc}j 



(3) 



The characteristic function (3) is computed like (2) applied to the real SaS 
random vector = . . . ,X],Xj) where Xj = X] + iX]. In this 

paper, when X"^ is complex we denote Txd the unique spectral measure corre- 
sponding to the real vector X"^*^ = {Xl, Xf, . . . , Xj, Xj). Let be the Fourier 
transform of the spectral measure r^d. When X'^ is real, it is given for all 

e = iei,...,ed)eR'' by: 



cos 



Sd)- 



(4) 



3 



For the complex case, it is given by, 

m= [ cos (nej20^s^^dT^^, (5) 



■2d 



where Sj = s] + is] and 9 = {9i, . . . ,9^) e 



Remarks 2.1 The finite measure T x<i is defined on the unit sphere which is 
compact this implies that it have finite moments of all orders, by the sequel 
is, in particular three times differentiable and we have: 



yid9jdak 



d9MM. ^ 



Od) = SiSjSk sin ^T^e 9eSi^ ^ dV ^ 



(6) 



9j is real then these operators are the usual partial derivatives. 

The covariation was conceived to replace the covariance which is undefined 
for SaS variables. It is given in the real case by: 

[x,,X2U= f si4"-i>rfr(x„x2)(5i,52), (7) 

where s^^'* = sign(s)|s|'^, see Samorodnitsky and Taqqu 19(]. Similarly for 
Xi = Xl + tXf and X2 = X^ + the covariation of Xi on X2 is given by: 



[Xi,X2], = / {s\ + tsl){s\ + isl)<^-^>dV^x, 

J§4 



X2){s\, si, s\, si) 



and z'^^^ = \z\^~^l, see Cambanis [l!]. The last author have shown that inte- 
grals in (7) and (8) remain unchanged if we replace T(^Xi,X2) by a higher order 
^{Xi,...,Xa)- This property is easily derived from the fact that, for all complexes 
(aj) and (bj) we have: 



d d 



+ is'))^''' ^dTxc^. (9) 



It was also shown that, for 1 < a < 2, the map X 1 — > \\X\\a = i[X,X]a)^ is 
a norm on the vector space induced by the set of SaS random variables. For 
more details and further properties see Cambanis 



2. 1 Sufficient condition for the additivity of the covariation 



In this section we are interested in the additivity of the covariation function 
with respect to its second variable, that is [F, aXj + bXj\a = [Y,aXi]a + 
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[Y, bXj]a, for all a and b. This property is usually provided by the independence 
of Xi and Xj or more generally by: F {(si, s^) G E>d,Si.Sj 7^ 0} = 0, see 
Cambanis and Miller [3|; The idea of this paper is to introduce a new 
condition on the vector X'^ more general than the condition of independence 
and permitting the additivity of the covariation. The next theorem provides 
an answer of this issue. 

Theorem 2.1 For the covariation to be additive with respect to its second 
variable that is, for all io G {1, . . . ,d}: 

V6'i, . . . ,9d G C, [Xjg, 6'iXi + ■ — h 9dXd]^ = [XjQ, 9iXi]^ H h [Xi^^,9dXd]^. 

(10) 

It is sufficient that, for all i,j and k G {1, . . . ,d} not all equals, the Fourier 
transform fulfill the next condition: 



Examples 

• Let us remark that the result (11) of theorem 2.1 is a generalization of the 
independence condition. Indeed, suppose that Xi, . . . , Xd are pairwise inde- 
pendent real SaS variables, then according to [Samorodnitsky and Taqqu 
ll9[, p. 68], the measure Txd is concentrated on the intersection points of 
the sphere E>d and the base axis of Mf^. We can then write as: 

= «l['^(i,o,....o) + ^(-1,0,... ,0)] + ■ ■ ■ + ad[\o,o 1) + '^(0,0,...,-l)]' 

where ai, . . . , are positifs weights. In this case, it is easy to see that the 
Fourier transform of F^d can be given by: 

d 

0(6'i, ...,9d) = ^aiCos{9i). 

i=l 

If we partially differentiate this formula two times we find a null function. 
This implies clearly that the function given in this example satisfies the 
condition (11). 

• Non trivial example: To get a simple example of SaS vector X'^ with spec- 
tral measure F^^ determined by its Fourier transform (p and verifying the 
condition (11), we choose of the form: 

^{9,,...,9d)=j2v^i^^)^ 

1=1 

where (/9j's are three times differentiable even real functions. By applying 
the Bochner's theorem, in order that be a Fourier transform of a finite 
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measure it is necessary and sufficient that it is positive definite. Then, for 
instance, one can take the v'ij's as three times different iable characteristic 
functions of random variables. 



2.2 Bimeasure construction and Morse Transue integral 



We recall that a stochastic process ^ = G M) is SaS if all finite subset 
of ^ is also an SaS vector. This is equivalent to the fact that all finite linear 
combinations of elements of ^ are also SaS random variables, see Cambanis 



jj, Samorodnitsky and Taqqu [19j. Let us consider the increments application, 
mapping each interval [s,t[ to the complex SaS random variable dC, defined 
by: 

da[s,t[)=^t-^s. (12) 

According to Cambanis [H, under the conditions (I) and (II), the application 
dC, may be extended to a complex random measure on the Borel cr-algebra 
i3(]R). For the construction of our bimeasure we will need the following condi- 
tion. 

Condition : 

We will say that dC, satisfy the condition A if and only if for all n > 2, 
and for all, pairwise distinct Borel sets Ai,...,An, the SaS random vector 
{d^{Ai), ...,d^{An)) satisfy the additivity condition (11) of theorem 2.1. 

Suppose now that the random measure d^ satisfy the condition A and consider 
the complex valued set function F defined on B{W) x i3(M) by: 

' ^ ' . (13) 

(A, 5) ^ [daA),daB)]^ 

The application F is additive with respect to its two variables: it is a bimea- 
sure. Indeed, the additivity of the ffist component comes from the linearity of 
the covariation with respect to its first variable see Samorodnitsky and Taqqu 



19l |. For the second component, let Bi and B2 be two disjoint Borel sets then 



for all fixed A we have, 

F{A, Bi U B2) = mA), d^{Bi U B2)]a = m^), d^{B,) + ^^(^2)],. (14) 
Since d^ satisfy the condition A then, 

F{A, fii U B2) = mA), daB,)U + mA), daB2)]a = F{A, + F{A, B2). 



^ The letter A is for "Additivity" 
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We can show that F is also cr-additive with respect to its two variables. Indeed, 



according to Rao [18|] or Chang and Rao j6||, it suffices to show that, for 



all Borel sets decreasing to 0, -F(y4„, A^) converges to 0. Since ^ is right 
continuous with respect to the weak topology of convergence in probabilitjo, 
we deduce that if (v4„)„ decreases to then d^{An) converges in probability 
to 0. This implies that F(y4„,y4„) = ||(i,^(y4„) ||° converges to 0. 



In the case of second order processes, it is known that the bimeasure associ- 
ated with the covariance function is positive definite. This property plays an 
important role in the construction of integrals with respect to a bimeasure see 
Chang and Rao Q . In our case of SaS processes, the bimeasure F defined in 
(13) have the following property: 

Proposition 1 The bimeasure F defined in (13) have a property similar to 
Bochner's positive definiteness: for all complex zi, . . . ,Zn and for all pairwise 
distinct Borel sets Ai, . . . , An, we have: 



EE 

i=i j=i 



Zi[Zj 



,<a-l> 



F{Ai,Aj) > 0. 



(15) 



The proof of this property is easy, it suffices to use the condition A. Indeed, 



i=l j=l 



^<-'>F{A,A,) = T.T.^^i^^)^''''ndaA^),daA,)V 
i=i j=i 

n n 

Y,zdaAi),Y.z4aA, 

li=l 



i=l 

a 



Y.^idi{A,) 



> 0. 



i=l 



In all the rest of the paper we suppose that the bimeasure F satisfy a condition 
similar to bounded Frechet variations but weaker than Vitali'o. It is given 



^ The topology induced by convergence in probability is equivalent to the covari- 
ation norm topology, see [Samorodnitsky and Taqqu [191], p. 95]. 



^ In this paper we will not be interested in some complicated considerations, about 
different bounded variations nor in their relationship with the existence of different 
type of integration theory with respect to a bimeasure F, discussed in Chang and 
Rao 0]. We just recall that, F have bounded Vitali variations if it verifies: 



sup < 



^^\F{Ai,Bj)\,Ai X Bj distincts > < oo. (16) 



=1 i=i 
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by: 



sup ^^Y.^ai{aj)<''-^>F{Ai,Aj),Ai e S(R) and \a,\ < l| < oo. (17) 

Let B he a. fixed Borel set and consider the complex measure defined by, 
Fb : A I — > Fb{A) = F{A,B). According to (17) tlie complex measure Fb 
have a finite total variation on B{W) (for the definition of total variation see 
[Dunford and Schwartz |9|], p. 97], we will recall it later in (20)). As in Rao 
18|, we define, in the sense of Dunford and Schwartz j^, the integral of a 



bounded function / with respect to Fb that we denote by: 

h{f,B)^ f f{X)FB{dX). (18) 

Now, let / be a fixed bounded map and consider the set function B i — > 
Ii{f,B) which is also a complex measure with bounded total variation on 
i3(]R). Similarly, we define the integral of a bounded function g with respect 
to/i(/,.): 

h{f,9)= I g{\)hU.d\). 



In the same way, by integrating first with respect to the complex measure 
Fa{B) = F{A,B) and A fixed, one can construct another integral l2{f,g) = 

/ f{X)l2{dX, g) where I2{A, g) = / g{X)FA{dX). Note that in general the 

Jr Jr 

two inte gra ls Ji and I2 are not equal (see counter example in Morse and 
Transue [l6|]). When Ii{f,g) = l2{f,g) the couple {f,g) is said Morse Transue 
integrable (MT-integrable) and the common value will be denoted as a double 
integral: 

I{f,9)= I I f{X)g{X')F{dX,dX'). (19) 



As it was mentioned by Rao [18|, it is not easy to identify the class of all MT- 
integrable functions. It should be noted that complexes measures on i3(]R^) 
are particular bimeasures, in this case MT-integration is the same as the usual 
Lebesgue (Radon) integration. For simplicity, we will be interested in a partic- 
ular class that we construct as follows: consider u : A 1 — > i^{A) = E(f (rf^. A)) 
where v{d^, A) is the total variation of the random measure dC,. It is defined 
in [Dunford and Schwartz \d\ , p. 97], for all Borel set A by: 

v{d^,A)= sup |^|de(^,)|,(A,)^^, partition of a|. (20) 
/ finite I iei J 

According to Dunford and Schwartz [9|, the total variation v{d^, .) is a positive 
random measure. We deduce then that z/ is a positive measure. We denote by 
Aa{dC,) the completion, with respect to the L^ii') norm, of complex simple 
functions with bounded supports. 
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3 Spectral representation of the covariation function 



Cambanis [jj have given a spectral representation of the covariation function 
of processes X given in (1) where ^ have independent increments. This last 
condition is the key of his spectral representation. Now that we have given a 
weaker condition of additivity in theorem 2.1, we generalize the Cambanis's 
integral representation of the covariation to processes X where increments of 
^ are not necessarily independents but verifying the additivity condition A. 
Our spectral representation is in the sense of the general Morse and Transue 



16| integration with respect to the bimeasure F. 



3.1 Integral representation of the covariation function 



It is easy to see that for simple functions / = ^ fi^Ai and g = ^ (7^11^^ from 

i=l j=l 



hif, g) = W, ^7) = E E f^9JF{A, A,) = J fd^, J gd^ 



i=i j=i 



(21) 



Thus, simple functions of compact supports are MT-integrables. In order to 
generalize the property (21) to all functions of Ka{dC) we will need the forth- 
coming proposition: 

Proposition 2 Suppose that d^ satisfy the additivity condition A, then we 
have the properties: 

(1) For all A G B{M) we have, \\d^{A)\\a < '^av{A) where \E'q, is a constant 
depending on a. It is given by = r(l — -)■ 

(2) Let B e B{M) he a fixed Borel set. If A is verifying viA) = then the total 
variation of the complex measure Fb in A is null, that is v{Fb,A) = 0. 
This result is also true for the measure Fa{B) and B fixed. 

(3) Let B G B{R) be a fixed Borel set, then for all bounded function f G 
Aa{d^) we have the inequality: 

I / fdFB\<^o.mBm-' f \f\di^- (22) 

(4) Let f G Aa{dC,) be a fixed bounded function and denote by G{B) = 
I{f,B). then {^{B) = implies that the total variation of G vanishes 
I.e. v{G,B) = 0). 

The importance of the last proposition comes from the fact that it permits 
to get an integral representation of the covariation function with respect to 
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the bimeasure F. It transforms the covariation of two stochastic integrals with 
respect to ^ to a double integral with respect to the bimeasure F which can 
be seen as a covariation of increments of ^. 

Proposition 3 Let f and g be two bounded functions in Aa(c?^). The covari- 
ation of the stochastic integrals / fd^ on / gd^ is given by: 

[ fd^J gdA =11 /(A)(^7(A'))<"-'>F(rfA,rfA'). (23) 

The term in the right hand side is Morse Transue integral with respect to the 
bimeasure F . 

The equality (23) is similar to the Cramer integral representation of the co- 
variance of second order processes, see for instance Chang and Rao [6|. This 
result plays a fundamental role in this work because it permits to characterize 
the law of the stochastic process X by the bimeasure F when the determinis- 
tic functions /(t, .) are known. This property will be used later to classify the 
dependence structure of harmonisable SaS processes according to the mor- 
phology of the bimeasure F. 

Proposition 4 Suppose that ^ is a real or complex isotropi^^ a-stable process 
where dC, satisfy the additivity condition A. Then the bimeasure F defined 
in (13) is the unique bimeasure characterizing X and verifying the integral 
representation (23). 



3.2 Applications to harmonisable processes 



In this subsection we are interested in harmonisable SaS processes. This means 
processes having the integral representation (1) where ( /(t, A) = e**'^ ) and d^ 
verifying the condition A. We also suppose that the bimeasure F given in (13) 
is of bounded variation in Vitali's sense. In this case, according to Horowitz 



ll| or Kluvanek [IJ], F can be extended to a complexe measure on the cr- 
algebra i3(M^) induced by BiWj x i3(M). By the sequel the MT- integration is 
the same as the usual Lebesgue (Radon) sense, see for instance Chang and Rao 
0. Let us consider the covariation function denoted by C : (s, t) \ — ^ [X^, Xt\a 
that we suppose uniformly continuous, then: 

• We will say that X is covariation stationary if the function C{s,t) 



^ A complex random variable Y is said isotropic if for all a;, the random variables 
Y and Ye"^ have the same distribution [Samorodnitsky and Taqqu [l9l |. p. 84 ]. It 
is known that isotropic random variables are parametric because their distribution 
is completely determined by their scale parameter (covariation norm). 



10 



depend only on the time difference s — t. 

• We will say that the process X is periodically covariated with a period 
T > if C(s + T, t + T) = C{s, t) for all s and t. 

• We will say that the process X is almost periodically covariated if C 
is an almost periodic function in the sense of Bohr. Dehay and Hurd ^] or 
Dehay 0. 

Remarks 3.1 It should he noted that with Cambanis's spectral representation 
where C have independent increments, the covariation of harmonisable pro- 

cesses .s e.-pressei as C(M) = / .here MA) = mWK = 

[d^{A), d^{A)]a is a spectral measure on B{M.). We see clearly that this kind 
of processes are always covariation stationary. In our general spectral repre- 
sentation, the covariation function and hence the dependence structure of the 
process depend on the structure of the bimeasure F. Indeed, by using (23) and 
the fact that ([e^^^ ^<a-i> _ g-«tA covariation function can be written as: 

i{s\~t\') 



C{s,t) = / / e'^'^''^^F{d\,d\'). 

jR JR 

The importance of our representation in the harmonisable case is that the 
covariation function is expressed as the Fourier transform of F. This result is 
very interesting: firstly it provides valuable informations on the link between 
dependence nature of the process and the morphology of its bimeasure. Secondly 
it makes possible statistical estimation of F from the process observations. 



According to Gladyshev [10|, the periodicity (resp. the almost-periodicity) 
of the bivariate covariation function C is equivalent to the periodicity (resp. 
the almost-periodicity) of the univariate maps, Cr '■ t i — > C{t -\- T,t) for all r 
fixed. As for periodically and almost periodically correlated processes see Hurd 



13l : Il2l |. this suggest that for all fixed r the function Cr may be expressed as 



a Fourier series. Indeed, when is periodic of period T, 



aw-E^^Me*"^*- (24) 

feez 

The convergence sense of the right hand series depends on the smoothness of 
C-r- The terms afe(r)'s are the order Fourier coefficient given by, afc(r) = 

— / Cr{t)e~^~^dt. Similarly, in the almost periodic case, the Fourier decom- 
T Jo 

position is given by: 

a(t) ~ E (25) 

I rM 

where = lim — / Cr{t)e~^^^dt are Bohr-Fourier coefficients of and 

^ M^oo M J-M 

Jt = {l & such that, a^(r) ^ 0} is at most countable, see for instance 
Dehay and Hurd In this last paper there was shown that according to 
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the uniform continuity of C, the Bohr- Fourier coefficients a^(r) are con- 
tinuous with respect to r. This fact permit to show that = \Jt >^r = 
{7 G M, such that, a^(r) 7^ for some r} is also countable, see Dehay and 
Hurd 

In the next proposition we show that some fundamental results concerning 
spectral representation of periodically and almost periodically correlated pro- 
cesses remain true for their analogues periodically and almost periodically 



covariated a-stables processes. As in the works of Hurd [12; ll3[, we give a 
necessary and sufficient condition for X to be periodically or almost periodi- 
cally covariated or covariation stationary. 

Proposition 5 (1) The process X is covariation stationary if and only if 
the bimeasure F is concentrated on the diagonal line. 

(2) The process X is periodically covariated with a period T if and only if 
the bimeasure F is concentrated on the lines parallel to the diagonal and 
equidistant of that is the lines {Lf^)^ defined by, = {(A, A') G M x 
M, /A-A' = ^}. 

(3) the process X is almost periodically covariated if and only if the bimeasure 
F is supported by lines parallel to the diagonal but not necessarily equidis- 
tants i. e. the lines {L^)^ where, L^ = {(A, A')gMx]R, /A — A' = 7G J} . 
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4 Appendix 



Proof of theorem 2.1: 

We begin with the case when X"^ is an SaS real vector and then we generahze 
to the complex case. 
Case when X'^ is real: 

We begin with demonstrating the additivity in the case of three components. 
By definition of the covariation and by using the result (9) it is easy to see 
that, for all io,j and k such that j ^ k and for all reals 9j and 9^ we have: 

= [Xia, 9jXj + 9kXk]^ - [Xi^^,9jXj]^ - [Xig, 9kXk]^, 

= ^ Sio {{9jSj + 9k.Sk)^" - {9jSj)^" - {9kSk)^" dVxd, (26) 

where A(sj, Sk) — {9jSj + ^feSfe)^" — {9jSj)^°' — (^feS^)^" For the rest 
of the proof we will need the next lemma: 

Lemma 6 Let 1 < p < 2, then for all real s, we have: 

' r(l-p)cos(f)yo tP ^'Jo tP 

where F is the usual gamma function. 
With the equality (27) we deduce: 

A(.„ ..) = sin((g,., + 9,s,)t) -smi9,s,t) - sin(g,.,t) ^^8) 

By using the classical trigonometric properties, 



sin(p) - sin(g) = 2 cos(^) sin(£^), 
cos(p) - cos(g) = -2sin(2±2) sin(2^). 



(29) 



it is easy to see that: 

Jo t"' 
We replace A{sj, Sk) in (26) and apply the Fubini theorem, which is applicable 
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in our case, we have: 

{ ...sm(%^*)sin(^*)sin(fe±^t)dr,. 
fl=-4.„^ '-^ ? ^- ^ *, (31) 

^ Jo 

where S(x,y) = [ Sj„ sin(— ^) sin(^-^) sin( — "1^ ^ '^ )drxd. Since the unit 

j§d 2 2 2 

sphere is compact, then the measure T^^ have finite moments of all or- 
ders. By the sequel the function S{x, y) is infinitely differentiable. Its partial 
derivative withe respect to y is given by: 

a^^^'^^ = 2 X/io«feSin(^)sin(^^-^)drx.. (32) 

By differentiating the above equality, with respect to x, we have: 

d'^S 1 f 

= 4 jg SioSjSfcSin((xSj +ysfc))drxd. (33) 

Let us remark by using (6) that: 



q2c 1 a'iA 

Therefore, according to the condition (11), we have: 

V^^^fcGM ^(^„^fc) = 0. (34) 

By (31) and (34), in order to show that, [Xi.OjXj + OkXk]^ = [Xi.OjXj]^ + 
[Xj, ^fcXfc]^, it is sufficient to show that (34) is equivalent to the fact that 5" is 
a vanishing function, that is: 

Vx, y, {x, y) = <J=^ yx,y,S {x, y) = 0. (35) 

The second imphcation is obvious because if 5" is a null function, it is the same 

for its second derivative. Now assume that we have, Vx, y, (x, y) = 0. We 

oxoy 

dS 

deduce then that \/x,y, ^—(x,y) does not depend on x. This implies that, 

oy 

dS dS 

\/x, y, -^{x, y) = -K~(0, y)- By replacing the right term of this last equality in 

dS 

the formula(32), we find that \/x, y, -^{x, y) — 0. In the same way, this implies 

y 

that S{x, y) does not depend on y, consequently Vx, y, S{x, y) = S{x, 0) and 
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by replacing this in the formula of S, given in (31), we see that S is null 
function. 



The proof of the general case is done with a technique similar to the one of 
three variables discussed above. Indeed, let 61,62, ... ,0d be real numbers and 
fix io e {l,..,d}. We note D the term expressed as: 



D 4 [X , + • • • + 6,X,l - [X„, ^^,XJ„ [X , ^,XJ„, 

= [X , ^^,X, + . ■ ■ + 6M. - [X.,,6,X, + ... + 6,_,X^J^ - [X , ^,XJ„, 

+ [X , ^,X, + • • • + 6,_,X^_,l - [X , ^,X, + • • • + 6,_,X^_,l - [X^^,6,_,XX, 

+ ■■■ + [X , e,X, + 6,X,l - [X , e,Xj„ - [X , e,Xj„, 

j=0 

(36) 

For all j e {0, . . . ,d — 2}, the terms D/s are given by: 

- [^zo'^i^i \- dd-j^d-j]c, - [X^.^^i^i H \- 6^_._j^X^_._^]^ - [Xi^,6d-jXd-j]„, 

= [ Si, [{6,si + --- + 6,_.s,_.y 



, <a-l> 



//I n ~.<a — l> , .<L-^ — 1>] 

-(^isi + --- + ^,_^_is,_^_j -[6^_.s^_.) 

(37) 

Like in the case of three variables discussed above, we use the equality (27) of 
lemma 6 and then we apply the Pubini theorem. We find an equality similar 
to (26). It is given by: 

K sin(g^t) sin(gi^i±::l%^i^t). sin( ^"^+-+;'----^^-' t)dr^d 

Dj = dt, 

_ f+o° A,{te,,...,t6a_,) ^^ 

(38) 

where. 



d 



, ^.-rVA • d.s, + --- + 6,_^_^s,_._, 
sm — sm — — 



2 



f6,s,-\ ^6, .s, .\ 

sin ^ '-' jdr^,. 



(39) 

First, the functions Aj's are two-differentiables because the measure Txd have 
finite moment of all orders. With the same technique as in (32) and (33), for 
all k — 1, . . . ,d — j — 1, by differentiating the function with respect to 6^ 
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and 9d~j-i we have: 



QQ^{e.. . . . , Od-,) = - sm( )dT^., 

(40) 

By using (6) and the condition (11), we deduce that: 

Vfc = 1, ■ ■ ■ , d - J - 1, V^i, ■ ■ ■ , 9a-j, {Oi, ■ ■ ■ , e^-j) = 0. (41) 

According to (38) and (41), in order to show that D = 0, it suffices to show 
the next equivalence: 

Vfe ^ 1, ■ ■ ■ , d - J - 1, V^i, ■ ■ ■ , da-j, (^1, ■ ■ ■ , 9a-j) ^ 

h d i 



V6'i, . . . , 9d-j, ^j{9i, . . . , 9ci-j) — 0. 



A trivial implication is that if Aj is a null function, it is the same for its 
second derivative. For the other implication, let us assume that \/k — 1, . . . , d— 

d'^A- 

j — 1 and V^i, . . . , 9d-j we have, ^ — (^i, . . . , 9d-j) — 0. One concludes, 

therefore, that V6'i, . . . , 9d-j-i the function (6^1, . . . , 9d-j) does not depend 

(j9d-j 

dA- 

on ^1, ... , 9d-j-i, which allows us to write, V^i, . . . , 9d-j, {9i, . . . , 9d-j) — 

o9d-j 

dA- 

^ (0, . . . , 0, 9 d-j)- By replacing the right term of this last equahty in the 
(j9d-j 

dA dA- 
formula ^ , we find that V6'i, . . . , 9d-j, ^ (^i, • • • , 9d-j) = 0. It also 
d9d-j o9d-j 

implies that does not depend on 9d-j, consequently V, 9i, . . . ,9d-j, we 
have the equality Aj{9i, . . . , 9d-j) = Aj(9i, . . . , 9d-j-i, 0). It is sufficient to 
replace 9d-j by in the formula (39), to see that is a vanishing function. 

Cage where X'^ is a complex vector: 

The prof of theorem 2.1 when X'^ is complex is similar to the real case. We 
give an idea on how to show this result in the case of three variables. First 
by definition of the covariation, for all io,j and k not all equals and for all 
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complex numbers 6j and 6k we have: 



= / ({OjSj + OkSk)^'' - (OjSj)^" - [OkSk)'^" dTx2d, (42) 

= / SioA(Sj-,Sfc)dr; 



X2d, 



where A{sj,Sk) = {OjSj + OkSk)^" — {djSj)^" — {dkSk)^" ■ Instead of 
using lemma 6, we use the following result: 

Lemma 7 Let z = a + ih he complex number and p positifs number such that 
< p < 2. Then by using the complex notation x = s + it, we have the 
following equality: 

roc roo I - cos(TZe(xz)) , , , M ,„ X 



-P i|i — cos(^ — j /■'^'^ £ 

where c(p) = 2~ — / 11 + sin(2^)| 2 d6. We have a result anal- 

p Jo 

ogous to (27): 



z 



= loo i-oo (a;)<^+^> ^^^^ = i-oo i-oo li)^^^^^'- 



sin(7?.e(xz)) , , a sin(7?.e(x2;)) 

(44) 

In this case, by applying the equality (44), for all x = Xi + 1x2 we have : 

/■ /■+°° sin(7^e(6'jSj + 9kSk)x) — sin(7?.e(6'jSjx)) — sin(7?.e(6'fcSfcx)) 
A(Sj,Sjj = Pa / 73<7^> 1 



(45) 

The rest of the proof is the same as in the case of real SaS vectors, where we 
use the lemma 7, the trigonometric properties (29) and the complex derivative 
operator as in (40) and (41). Finally the same reasoning as in the real case 
shows the results. 

Proof of proposition 2 : 

(1) Let A G B(R) be a fixed Borel set. Since d^{A) is a real (resp. isotropic 
complex) SaS random variable then according to Cambanis and Mi- 
amee ^ or Samorodnitsky and Taqqu [l^ , the p^^-order fractional mo- 
ments of d^{A) are given by E\d^{A)\P = Sa{p)\\d^{A)Z where S^ip) = 

2^ r(i_£)r(i) (^^SP- 'S'q(p) = r(i-g) ° )• In particular, for p = 1, 
Sail) < ^1 then |M^(A)|U < ^„.E|rf^(A)|. According to the definition 
of the total variation, defined in (20), it is easy to see that |(i^(y4)| < 
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v{d^,A). Therefore, we have: 

\mA)\\^ = ^„E|(ie(A)| < ^^E{v{d^,A)) = ^MA). (46) 



(2) Let B G B{M.) be fixed Borel set, then by definition of the total variation 
of the complex measure Fb given in (20) (see Dunford and Schwartz [9^) 
and replacing Fb by its expression we have: 

v{Fb, A) = sup < ^ |Fb(A)|, {Ai)i(zi partition of A L 
/ finite I ie/ J 

= sup 1^ |[(i^(A),(i^(5)]a|,(^i)ie/ partition of aI. 
/ finite I ie/ J 

It is known that \[d^{Ai), d^{B)]a\ <i[c?^(A)IUIM^(A)||a"S see for in- 
stance [Samorodnitsky and Taqqu [l9|, p. 96]. Then by using (46) we 
deduce that, 



v{Fb,A) <\\daB)ra-' sup ^ ||rfe(^.)IU,(^.)ie/ partition of A , 

/ finite lie/ J 

< ^,||rfe(5)||r' sup \Y.u{Ai), {Ai),^j partition of a| , 
/ finite I ie/ J 

<^jdm\\T'^iA). 

By the sequel, if i^{A) is null then v{Fb, A) = 0. 
(3) First we begin by showing the inequality (22) for simple functions. Indeed, 

n 

let / be a complex simple function defined as / = ^ fi^Ai, then : 



i=l 



fdF^ 



< mm 



j2fimA),daB)i 

i 



Y^fdi{Ad,di{B) 



< v[/J|rfe(5)|ir El/^l-E(^(rfe,^.)), using (46) 



a-l 



\f\du. 



(47) 

If now / is a bounded function in Aq,(^), then there exist as sequence 
of simple functions (/„) such that J \ fn — fldu converge to as n tends 
toward infinity. First let us remark that the sequence {J^ fndFB)n is con- 
vergent because it verifies the Cauchy condition. Indeed, according to the 
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inequality (47), for all integers m and n, 



/ fndFB - I fmdFs <^a\mB)\C I \fn-fm\dv. (48) 

On the other hand, since {fn)n converges to / in Liiy) we can then 
extract a subsequence {fn^)k that converges to / (//-almost everywhere). 
This means that there exists a Borel set D G i?(]R) such that z/(R\D) = 
and fuk^D converges toward /.Id as k tends to infinity. Using the 
second result of these proposition and the fact that viM. \ D) = we 
deduce that f (Fe,M \ D) = 0, this implies that the subsequence {fnjk 
converges v{Fb, .)-almost everywhere to /. The functions /„j. and the 
measure v{Fb, .) are bounded, then by applying the Lebesgue convergence 
theorem in the case of Dunford-Schwartz integrals [Dunford and Schwartz 
[sl, p. 151], we deduce then that the subsequence {J^ fnkdFB)k converge to 
J^fdFs- Using the uniqueness of the hmit we deduce that the sequence 
U^fndFB)n converge to J^fdFs- Finally, by using (47) we have: 

I / fndFBl < \mB)\\T' I \fn\dy. (49) 

JR JR 

The inequality (22) is then obtained by passage to the limit in (49). 
(4) Now let / G Ka{dC} be a fixed function. Using the definition of the total 
variation of G, for all B G i3(]R), we have : 

v{G,B)= sup Ij] I G(5,)|,(5,),e, partition of b|, 
I finite I ie/ J 

= sup i /c/FbJ, (5i)jg/ partition of B I , 

/ finite I ie/ "^'^ J 

< / \f\dv sup \\di{B,)\\l-\ {B,),^i partition of b| . 

•'^ I finite I iei ) 

(50) 

By using the inequality (46) in (50) we deduce that: 

viG, B) < i^^r [ \f\di^- sup \J2 HBi)r'\ {B,\ei partition of b| . 

•'^ I finite Ue/ J 

(51) 

If i'{B) = then for all finite index set I and for all partition {Bi)i^i 
of S, the BiS are subsets of B, therefore v{Bi) = for all i and we 

have sup <^ WiBi)\''-\ {Bi)i^i partition of B [ = which impl 

/ finite I i&i ) 
that v{G,B) = 0. 

Proof of proposition 3: 

Let / and g be two bounded functions of Aa{dC,) and /« = X! fl^^^" 9- 



les 



20 



g^^Bv- two sequences of simple functions in Aa{d^) such that |/„ — fldu 

j 

and Julgn — gldiy converge to 0. The equahty (23) is satisfied for the sequences 
of simple functions {fn)n and {gn)n- Indeed, let m and n be two integers, then: 



/ fndC, / gmdC 



The Borel sets (-B")^ are pairwise disjoints and since d^ satisfy the condition 
{A) we deduce that: 

/ fnd^, f gmdA = EE/r fei^""'^ mA^),daBj')u 

-T.T.fr{9TT^''^F{A7,B^), (52) 

i j 

-II fn{X){gm{X')f F{dX,dX'). 

The result (52) remain true for any function / and g in Aa{d^). First, remark 
that: 



<a-l> 



EE/r fo"') [dm),daBT)]a = E (aT) / fn{X)FB-{dX). 

i j 3 

Firstly, for a fixed m and using the same reasoning as in the proof of the 
third assertion of the proposition 2, we have: j fn{X)FB-^{dX) converge to 

j f {X) Fbj^ (dX) . On the other hand, using (46) we can derive \\J^ fn — fd^\\^ < 
/ l/n ~ fldj^- We have then the inequality. 



fndi, / gmdi - / fdi, / g^d^ 

JR Ja LJM JR 



< 

< ^ 



fn- fd^, / gmd^ 

JR JR 

I fn-fdi 

JR 

gmdi 



a-1 



I gmdC 

JR 

|/„ - f\diy. 



Since (/„) converges to / in i^^(i^), we deduce that / /„rf^, / gmd^ 

UR JR 



con- 



verges to 



By the sequel as n tends to infinity in the equality 



fd^, / gmdC 

JR 

(52) we succeed to: 

/ fd^, / gmdC = E (9 

jR JR ja A 

= J {grr.{X'))<'^-'>G{dX'), 



'"^^""'^ ' f{X)FBrr{dX), 



(53) 
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where the measure G is defined in the proposition 2. Since {gm)m converges to 
g in Li{u). there exists a subsequence {gmu)k of (fi'm) that converges to g {v- 
almost everywhere). By applying the fourth assertion of proposition 2, we have, 
(fi'm(A'))^"~'^^ converges to (5'(A'))<"~^> (^(C; O'^^^iost everywhere). Finally, 
since G is of bounded total variations, then the Lebesgue convergence theorem 

ensure the convergence of j (5^(A'))<"~^>G(rfA') to j {g{X'))<''-^>G{dX').By 

a similar technique as for (/„) , 



fd^, / gmd^ 



fd^, / gd^ 



< 



fdC 



gm - gd^ 
fdi 



a-l 



gm - g\dy 



a-l 



Since (fc) couwrges m L'(.) to g then the covariafon, [/ / fa<i« 



converges to 
(53) we obtain. 



/ fdi, j gdC 



. Finally, as m tends to infinity in the equality 



/ fdC, [ gdc] = / {g{X'))<''-'>Gf{dX')=h{f,g). 

.JR JR J a J 



(54) 



fdC, / gd^ 

JR 



The same reasoning on the integral I2 permits to show that 

h{f,g)- In conclusion we have shown that the couple {f,g) is MT-integrable 
and we have (23). 

Proof of proposition 4: 

Let A and B be any two Borel sets in M and suppose that there exist an other 
bimeasure F' verifying (23). In this case, by choosing the particular functions 
f — 1a and = 1^ we have : 

mA),d^{B)]^ = 1 1 lA{X)lB{X')F{dX,dX') = 1 1 lA{X)lB{X')F'{dX,dX'). 

which imphes that ^(^4, B) = F'{A, B) and the uniqueness of F is proven. 
Now let, ti,t2, ■ ■ ■ ,tn finite sequence of reals. Let us show that the character- 
istic function of the SaS random vector, (X^^, . . . , Xt^) is characterized by F. 
For this purpose we distinguish the real and the complex cases. 

The real case: Using the definition of the characteristic function of real SaS 
random vector, we have: 



*eiXt^+-+e„Xt„(l), 
e^p{-\\eiXt, + --- + 9nXtXJ, 
exp{-[0,X,, + --- + 0^X,^,0,X,,+ 
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Replacing Xt by its expression given in (1), it is easy to see that, 



i=l 



1=1 



5:^, / /(ti,A)de(A),E^^ / /(^i,A)rfe(A) 

1=1 1=1 

n . n 



Since the functions f(ti,.) are bounded in Aci((i05 which is the same as 

n 

Y^9if(ti, .), then according to (23) of proposition 3 we have: 



i=l 



Li=l 



i=l 



„ „ n n 
i=l i=l 



This imphes that, 



$x.„...,x.„(^i,...,^n)=exp|- 1 1 (X:^./(t.,A))(X:^./(t.,A'))<"-^>F(rfA,rfA0| 

(55) 

It is clear from this last expression that the characteristic function is expressed 
directly by the known deterministic functions f{t, .) and the unique bimeasure 
F. 

The isotropic complex case: It is known that isotropic complexes a- 
stables random variables are parametric. Their characteristic function is en- 
tirely determined through the covariation norm see for instance Cam- 
banis Uchaikin and Zolotarev 20 1. Therefore, using the definition of the 
characteristic function of the random vector (X^^, . . . , X^^) we have: 

^Xt^,...,Xt„idl, ■ ■ ■ ,9n) = $6»iXt-^+ - +e„Xi„(l), 

= exp {-Co \\9iXt^ H h 9nXt„ ||°} , 

= exp {-Co [e^Xt, + ■■■ + OnXt^.e^Xt, + ■■■ + dnXtX) , 

where Cq is given by cq = — / | cos(0)|°(i0, see [Samorodnitsky and Taqqu 

I — I 27r JO 

191], p. 86]. Similarly as in the real case we have: 

{» » n n 
-Co / / (E^^/(^«,A))(^^,/(t„A'))<"-^>F(rfA,rfA') 
i=i i=i 



Xt^,...,Xt„\.^l^ 



Which achieve the proof of the proposition. 

Proof of proposition 5 : The proof of (2) of this proposition is the same 
as in Hurd proposition 1] for the case of periodically correlated processes. 
The proof of (3) is the same as in the almost periodically case see Hurd 
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13l . proposition 6]. In order to give the idea of the proof, we recall here the 
proof in the periodically covariated case. First, for all fixed r, by using the 
representation (23) and the fact that (e**'^ ^<a-i> _ |g?iA ^a-2^~itx _ ^-itx ^ 

have: 



On the other hand, according to the periodicity of the covariation function, 
for all > we have : 



1 ^ 



N 

1 ^JrJr 



1 



2N + 



I I e-^e^*(^-^')D^(A,A')F(dA,dAO 



1 

with, Dn{\ a') = ^ e'^^^^-^'\ It is the Fejer's kernel (see Kurd (l3|). 



It is easy to see that: 



N 



1 si A - A' = ^, 



- ^ 1 sin((iV+l)T(A-AO) . , w , 2,^ ^^^^ 
2Ar+l sm(T(A-A')) ^ 7^ r • 

As tends to infinity in (56), Diy{X, A') converge to 11^ where L = Uki^zLk = 

{2nk \ 
(A, A'), A — A' = — - >. Since F is of bounded Vitali variations, then by 
-'J 

the Lebesgue convergence theorem which is true in our case, we have for all t 
and r: 

" e*(*+^)\-**^'F(c/A,c/A') = / / e<*+")^e-**^']li(A,A')F(c/A,c/A') 



Applying the inverse Fourier transform we see immediately that F is concen- 
trated on L = Ufcez ^k- 

Reciprocally, by using the fact that F is concentrated on L we have: 
[Xt+r+T,X,^r]a= I I e^(*+^)\-*^'e^^(^-^V(rfA,rfA'), 

JR JR 

e'^'^^^^e-''^' e''^^^"^'HL{\ \')F{d\, d\'). 



However, for all (A, A') in S, T(A - A') belongs to 2nZ then e'^^^"^') = 1 by 
the sequel: 

[Xt+^+T,Xt+T]a= [ [ e<*+^)"e-**^'F(rfA,rfAO = [Xi+.,Xi]„, 
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which achieve the proof of the proposition. 



Proof of lemma 6: We have the following classic equality that one can 
find, for example, in Nikias and Shao [l7t]: 



1 - cos(st) 



dt 



r(l-a)cos(^f^ 



a 



(57) 



Since 1 < a < 2, for all real s, the map t 



sin(st) 



^ is integrable on [0,+oo[ 



because when t tends toward infinity we have the asymptotic approximation 



I sin(t) I 



and when t tends to 0, 



I sin{t) I 



We can then differentiate 



the equality (57) under the integral sign with respect to s. We have therefore 
the following equality : 



+ 00 



sin(st) 



dt 



,<a-l> 



an. 



r(l — a) cos( — ). 



which achieves the proof of this lemma. □ 

Proof of lemma 7: By using the complex notations z = a + tb et x = s + tt 
it is clear that lZe[xz) = as + bt and |x| = ^/s"^ + t"^. Therefore, the calculation 
of the integral (43) is equivalent to, 

roo /-oo X _ cos(as + bt) 



(^2 _|_ t'^yt 



dsdt. 



(5J 



To calculate the integral 1, for all a and b, we use the change of variables: 

1 




t 
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[{a + b)x + (a - b)y], 
■[{b - a)x + (a + b)y\, 



(59) 



a2 + 62 

First, the Jacobian of the transformation (59) is given by: 

a + b a — b 



J 



b — a a -\-b 
+ 6^ + 62 



+ 62' 



A simple calculus shows that: as + bt = x + y et s^ + t^ = ^^^p-^, this implies: 

1 — cos(a; + y) 2 

—dxdy, 

(60) 



oo roo 



2—{a^ + h- 



oo /-oo X _ cos(x + y) 



°° (a;2 + y'^)^^^ ^ 



dxdy. 
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By a change of variables to the polar coordinates and thereafter using (57), 
we have: 

°° l-cos(x + y) , , /•o<^ 1 - cos(r(cos(^) + sin(^))) , ,^ 
^+^dxdy = I I -j^, rdrde, 



— oo J— oo -|- y^) 2 Jo Jo 



r(l-p)cOs(f) . 2^ 



r(i-p^cos(f) /-s- 



(P+2) 

cos(^) +sin(^)|^^d^, 
r" 11 + sin (2^)1 2 d^. 

^0 



p Jo 

By replacing this last equality in (60), we obtain (43). 

Let us consider the derivative operator in space of complexes numbers C de- 

d ( d d\ 

fined for 2; = a + by, — = 7- 1— \. As in lemma 6 we differentiate 

oz \oa 00 J 

with respect to z the equality (43) under the integral sign, we find: 

roo foo Q / 1 - C0S(71C(XZ))\ , , , , , , , 

LLt.[ i.|U> " }'^''t=c(p)^j.r. (61) 

g 

A simple calculus of the derivative show that, t^UI^ = pz^^~^^ by using the 

oz 

d 

same technique we have, — cos{TZe{xz)) — —xsm{TZe{xz)). We then use this 
results in (61) we conclude that: 

/■oo ^oo sin(7?.e(a;^)) , , , ^ ^„ /^^x 
/-J-coli)^"='"W^--". (62) 

This achieve the proof of lemma 7. □ 
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